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Abstract 

We study nucleon-nucleon scattering on the lattice at next-to-leading order in chiral ef- 
fective field theory. We determine phase shifts and mixing angles from the properties of 
two-nucleon standing waves induced by a hard spherical wall in the center-of-mass frame. 
At fixed lattice spacing we test model independence of the low-energy effective theory by 
computing next-to-leading-order corrections for two different leading-order lattice actions. 
The first leading-order action includes instantaneous one-pion exchange and same-site con- 
tact interactions. The second leading-order action includes instantaneous one-pion exchange 
and Gaussian-smeared interactions. We find that in each case the results at next-to-leading 
order are accurate up to corrections expected at higher order. 
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I. INTRODUCTION 



There have been several recent studies on the su 



)ject of Jattice simulations for low-en ergy 

16]. 
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nuclear physics using effective interactions [l, 2], [j, 4, [5, If], Q, fj, |jj 
Here we present a study of nucleon-nucleon scattering on the lattice at next-to-leading order 
in chiral effective field theory for momenta less than or equal to the pion mass. This 
analysis is part of a long-term effort to put the formalism of chiral effective field theory on 
the lattice and use the lattice action to study few- and many-nucleon systems. This is the 
first of a pair of papers on lattice chiral effective field theory at next-to-leading order. In 
this first paper we calculate lattice phase shifts and mixing angles and fit unknown operator 
coefficients at next-to-leading order using the spherical wall method introduced in [171 ] ■ In 
the companion paper we use the resulting lattice action to perform Monte Carlo simulations 
of dilute neutron matter. 

The organization of this paper is as follows. We begin with a review of the effective 
potential for chiral effective field theory at next-to-leading order and simplifications which 
can be made at low cutoff momentum. We then discuss a general test of model independence 
at fixed lattice spacing. We illustrate this model independence explicitly by repeating all 
calculations using two different lattice actions. The first action uses instantaneous one-pion 
exchange and same-site contact interactions at leading order. The second leading-order 
action includes instantaneous one-pion exchange and Gaussian-smeared interactions. Both 



of these lattice actions were first introduced in 



131 ] . We compute next-to- leading-order 



corrections for these actions and find results accurate up to omitted higher-order terms. 



II. CHIRAL EFFECTIVE FIELD THEORY AT NEXT-TO-LEADING ORDER 



A. Effective potential 



In the following q denotes the t-channel momentum transfer for nucleon-nucleon scat- 
tering while k is the w-channel exchanged momentum transfer. At leading order (LO) in 

n □ 

the Weinberg power-counting scheme [181 . 1191 ] the NN effective potential consists of two 
independent contact terms and instantaneous one-pion exchange (OPEP), 
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The vector arrow in a signifies the three-vector index for spin. The boldface for r signifies 
the three-vector index for isospin. We take for our physical constants m = 938.92 MeV 
as the nucleon mass, m n = 138.08 MeV as the pion mass, f n = 93 MeV as the pion decay 
constant, and qa = 1-26 as the nucleon axial charge. 

At next-to-leading order (NLO) the effective potential introduces seven independent con- 
tact terms carrying two powers of momentum, corrections to the two LO contact terms, 



and instantaneous two-pion exchange (TPEP) [20 
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conventions of 
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24j. Using the notational 



24j we have 
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NLO 



V LO + AV^ + V^ + V^ p . 
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The contact interactions are given by 



AV (0) = AC S + AC T (ai • a 2 ) 



(5) 



V® = C lQ 2 + C 2 k 2 + (C 3 q 2 + C A k 2 ) (<?! • a 2 ) + iCr}- (at + a 2 ) ■ (qx £ 

+ C % [a x ■ q) (a 2 ■ q) + C 7 U ■ k) (a 2 ■ k) , (6) 



and the NLO two-pion exchange potential is 
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L (q) [(<? • ^i) (q ■ #2) - q 2 (^i ■ £2)] 
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ti ^ 1 I a 2 ; 2 1 V 4m ^ + q 2 + q 

L(q) = -^ml + qnn ^ mi + q2 _ q 
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B. Simplified form at low cutoff 

In this paper we consider low-energy nucleon-nucleon scattering at momenta less than or 
equal to the pion mass, m n . On the lattice the ultaviolet cutoff momentum, A, equals 7r di- 
vided by the lattice spacing, a. If we focus narrowly on calculations of two-nucleon scattering 
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on the lattice we can take any sufficiently small lattice spacing satisfying A ^> m n . However 
serious numerical difficulties appear at large A in Monte Carlo simulations of few- and many- 
nucleon systems. In some attractive channels we might find spurious deeply-bound states 
at large A. In other channels we must deal with short-range hard-core repulsion becoming 
more prominent and producing fluctuating sign or complex phase cancellations. The gen- 
eral connection between sign/phase oscillations and repulsive interactions has been discussed 



in the literature in several different contexts 
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291 ] . The severity of the sign/phase 



problem scales exponentially with system size and strength of the repulsive interaction. 

In order to avoid these difficulties we set the cutoff momentum A as low as possible for 
describing physical momenta up to m^. We take A = 314 MeV ~ 2.3m 7T , corresponding 
with a -1 = 100 MeV. This coarse lattice approach is similar in motivation to the continuum 
low-momentum potential V\ ow k derived using the renormalization group 3(J, |3lJ. There has 



been some discussion in the literature about the consistency of the Weinberg power counting 
scheme at high momentum cutoff, starting with the work of 32|, |33| and more recently 



34 



35 



36 



37 



381 ] . Due to the computational reasons explained above, we are prevented 



from reaching high momentum cutoff scales in our lattice calculations where alternative 
power counting schemes may be useful. 

For nearly all \q\ < A we can expand the two-pion exchange potential in powers of 
q 2 /{4ml), 



3 4m 2 
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This expansion fails to converge only for values of q near the cutoff scale A w 2.3m 7r , where 
the effective theory is already problematic due to large cutoff effects of size O (q 2 /A 2 ). From 
a practical viewpoint there is no need to keep the full non-local structure of V^^p at this 
lattice spacing. Instead we simply use 



V LO = V® + y° PEP , 



(12) 
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V NLO = V LO + AV^ + V^\ (13) 

where the terms in Eq. (TTTT) with up to two powers of q are absorbed as a redefinition of 
the coefficients AV {0) and V {2) . This same approach can be applied to the two-pion ex- 
change potential at next-to-next-to-leading order (NNLO) and higher-order ra-pion exchange 
potentials. 

C. Model independence at fixed lattice spacing 

The usual test of model independence in low-energy effective field theory calculations 
is to check for sensitivity on the cutoff scale A. The difference between calculations at a 
given order for two different cutoff scales K\ and A 2 should be no larger than the omitted 
corrections at the next order. On the lattice this test is problematic since we cannot change 
the lattice spacing by a large amount due to computational constraints. Fortunately cutoff 
independence is just one of many ways to test model independence. In the following we 
discuss a general class of methods to check model independence at fixed lattice spacing. 

Let us use the notation yQ n l kn to denote two-nucleon operators with the following prop- 
erties. yQ n / A " i s a sum Q f local two-nucleon interactions each with at least n powers of 
momenta and an analytic function of momenta below the cutoff scale A. We use the adjective 
"quasi-local" to describe yQ n l kn since the interactions are purely short range. Quasi-local 
operators of this type arise naturally in improved lattice actions of the type considered in 
13]. 

At fixed lattice spacing we consider two modified lowest-order actions of the form 

VL Ol = \/ 1 (0) + ^ OPEP + n Q2/A2 , (14) 
V L02 = V^ + V OPEP + V 2 Q2/A2 , (15) 

where V® ^ A and V® ^ are different quasi-local operators with at least two powers of 
momenta. The leading-order interactions are iterated nonperturbatively and so, when 
fitted to physical low-energy data, the contact terms and in general have different 
coefficients. However we expect low-energy physical observables such as scattering phase 
shifts, binding energies, etc., should agree up to differences comparable to the omitted NLO 
corrections. 
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Similarly at NLO we may consider modified actions of the form 

Vj^ = V Uh +AV® + vl» + V?' A \ (16) 
Vv L02 = V L02 +AV<V + V^ + Vp A4 , (17) 

where V^ 1 ^ and vf 1 ^ are different quasi-local operators with at least four powers of 
momenta. Once again low-energy physical observables should agree up to differences com- 
parable to the omitted corrections at the next order. 

At any given order the effect of small changes of the lattice spacing can be reinterpreted 
at fixed lattice spacing as a renormalization group transformation on higher-order local 
operator coefficients. In principle this type of modification is covered by our condition of 
model independence for general V^"^ A " . Nevertheless it is good to check explicitly as many 
variations as possible, and studies at different lattice spacings as well as different functional 
forms for yQ n /^ n are planned for the future. 



III. LATTICE FORMALISM 



A. Lattice notation and cubic symmetries 

In this paper we assume exact isospin symmetry and neglect electromagnetic interactions. 
We use n to represent integer-valued lattice vectors on a three-dimensional spatial lattice 
and either p, q, or k to represent integer- valued momentum lattice vectors. I = i, 2, 3 
are unit lattice vectors in the spatial directions, a is the spatial lattice spacing, and L is 
the length of the cubic spatial lattice in each direction. We use the Euclidean transfer 
matrix formalism defined in [3] with lattice time step a t , and the integer n t labels the time 
steps. We define a t as the ratio between lattice spacings, a t = a t /a. Throughout we use 
dimensionless parameters and operators, which correspond with physical values multiplied 
by the appropriate power of a. Final results are presented in physical units with the 
corresponding unit stated explicitly. As in [l^ we use the spatial lattice spacing a = (100 
MeV) -1 and temporal lattice spacing a t = (70 MeV) -1 . 

We use a and to denote annihilation and creation operators. To avoid confusion we 
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TABLE I: Irreducible S0(3, Z) representations 



Representation 


Jz 


Example 


At 


mod 4 


^0,0 


T X 


0,1, 3 mod 4 


{y 1)0 ,Fi,i,yi,_i} 


E 


0, 2 mod 4 




T 2 


1,2, 3 mod 4 




A 2 


2 mod 4 


Y 3 ,2-Y 3 ,-2 
V2 



make explicit in our lattice notation all spin and isospin indices using 

a 0,0 — a T,P> a o,i = a T,n> (18) 

ai,o = aj.,p, ai.i = a;,n- (19) 

The first subscript is for spin and the second subscript is for isospin. We use tj with 
I = 1,2,3 to represent Pauli matrices acting in isospin space and as with S = 1,2, 3 to 
represent Pauli matrices acting in spin space. We also use the letters S and / to denote 
the total spin and total isospin for the two-nucleon system. The intended meaning in each 
case should be clear from the context. 

On the lattice the rotational symmetry is reduced to the cubic subgroup SO(3,Z) of 
SO (3) while isospin symmetry remains intact as the full SU(2) symmetry. There are five 
irreducible representations of the cubic rotational group. These are usually written as Ai, 
Ti, E, T2, and A 2 . Some of their properties and examples in terms of spherical harmonics 
Yl,l z {9,4>) are listed in Table HI The 2 J + 1 elements of the total angular momentum 
J representation of SO (3) break up into smaller pieces consisting of the five irreducible 
representations. 

We use the eight vertices of a unit cube on the lattice to define spatial derivatives. For 
each spatial direction I = 1, 2, 3 and any lattice function f(n), let 

A i /(n) = i Yl (-ir + 7(™ + ^), *=i/ii + ^2 + i/ 3 3. (20) 

^1,^2,^3=0,1 

We also define the double spatial derivative along direction I, 

V 2 if{n) = f(n + 1) + f(n - 1) - 2f(n). (21) 
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B. Densities and current densities 



We define the local density, 



P at ' a (n) = Yl (22) 



which is invariant under Wigner's SU(4) symmetry 39J. Similarly we define the local spin 
density for 5 = 1,2, 3, 

P a s' a (n)= E 4/™) Ws^ai'jin), (23) 

ij,i'=0,l 

isospin density for / = 1, 2, 3, 

pf' a (n)= Y «L-(^)N,,,a^(n), (24) 

M,j'=0,l 

and spin-isospin density for S, I = 1, 2, 3, 

Psf(n) = E a h( n ) \. a s]u' [ T i]jf a i',j'(^)- (25) 

,j'=0,l 

For each static density we also have an associated current density. Similar to the defini- 
tion of the lattice derivative A; in Eq. ([217]) . we use the eight vertices of a unit cube, 

V= v x \ + u 2 2 + u 3 3, (26) 

for v\, i>2, z/3 = 0, 1. Let for I = 1, 2, 3 be the result of reflecting the / th -component of 

V about the center of the cube, 

u{-l) = u+{l-2ui)l (27) 

Omitting factors of i and 1/m, we can write the Z th -component of the SU(4)-invariant current 
density as 

Tlf> a (n) = \ J2 E + *H))«u(* + V). (28) 

v\, V2, f3=0,l ij=0,l 

Similarly the Z th -component of spin current density is 

= \ E E (-^H^ + *(-')) Nii' + (29) 

^1,^2,^3=0,1 jj,i'=0,l 

Z th -component of isospin current density is 

= | E E (-^ l+ %(n+u(-l))[r I ] jjl a i ,Afi + P), (30) 

Vl,V2, ^3=0,1 i,j j'=0,l 
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and / th -component of spin-isospin current density is 

n &» = 7 E E (- i r i+ %^ + ^-n)Ws] i Ari] jr ^An + ^. (31) 



^1,^2,^3=0,1 i,j,i',j'=0,l 



C. Instantaneous free pion action 



The lattice action for free pions with purely instantaneous propagation is 

S^im) = + 3) ^ 7Tj(n, n t )7Tj(fi, n t ) - a t 7r/(n, n t )7r 7 (n + f,n t ), (32) 

ii,nt,I n,nt,I,l 

where 7T/ is the pion field labelled with isospin index I. We note that pion fields at different 
time steps n t and n' t are decoupled due to the omission of time derivatives. This generates 
instantaneous propagation at each time step when computing one-pion exchange diagrams. 
It is convenient to define a rescaled pion field, rc'j, 



q w = a t (ml + 6). 



Then 



SVttW) = \ E ^ii.n,n t )n' I (n,n t ) - — E ^1(^,12^^(71 + 1, 



n tj 



n,nt,I 

In momentum space the action is 



n,nt,I,i 



i,k 



1 _ Ut 

2 q n 



(33) 
(34) 

(35) 
(36) 



The instantaneous pion correlation function at spatial separation n is 

//^ \ / /j-? A J Z^7Tj 7i'j(n, n^TTjiO, n t ) exp [— 5' 7r7r ] 
717(71, ritj-Kj^, n t ) ) = c f — - — ^ (no sum on I) 



L 3 



/ Dn'j exp 

e- i2 ^ n D n (k), 



where 



l-^E.cos(^)- 
It is useful also to define the two-derivative pion correlator, Gs 1 s 2 {n)i 

Gs^in) = ( A 5l 7r^(n,n t )As 2 717(0, n t )\ (no sum on /) 



(37) 



(38) 



Tq E E (-ir Sl (-l)^(<(^ + ^-^n t K(0,n,)). (39) 



^1,^2,^3=0,1 v ' x ,1/^,1/^=0,1 
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IV. LATTICE TRANSFER MATRICES 



A. Leading-order actions LOi and LO2 



The analysis in 13] considers the lattice path integral with and without auxiliary fields 
as well as the lattice transfer matrix with and without auxiliary fields. All four formulations 
are shown to be exactly the same. Here we discuss only the lattice transfer matrix without 
auxiliary fields. This formulation is the most useful for calculating nucleon-nucleon scat- 
tering phase shifts and mixing angles. In simple terms the Euclidean transfer matrix is the 
exponential of the Hamiltonian exp(— HAt), where At equals one temporal lattice spacing. 
For example we can write the free-nucleon transfer matrix as 



M free =: exp (-H {rce a t ) : 



(40) 



where the :: symbols indicate normal ordering. We use the 0( a 4 ) -improved free lattice 
Hamiltonian, 



free 



49 
12m 



Y Y a lj(n)a itj (n) 



ft i,j=0,l 



—Y y y 



ft i,j=0,l Z=l,2,3 



alj(n)aij(n + 1) + alj(n)a ii j(n - I) 



40m 
1 

180m 



Y Y Y rlji^aijin + 21) + a\ d (n)a id (n - 21) 



ft i,j =0,1 1=1,2,3 



(41) 



ft i,j =0,1 1=1,2,3 

With the interactions included we take the standard lattice transfer matrix defined in 
131 for LOi, 



M LOl =: exp < -H {rcc a t - ^Ca t y Jp ot '°(fi) - -Cpton y y |pf ,a 



1 2 



n 



1 ft 



+ 



8/X 



y y G Sl s 2 (ni - "2)p^' i l(ni)p5; ) j(n 2 ) 



(42) 



Si,S2,I n\,n2 

where C is the coefficient of the Wigner SU(4)-invariant contact interaction and Gp. is the 
coefficient of the isospin-dependent contact interaction. For C and Cp we use the values 



C = (3C I=1 + C I=0 ) /A, 
10 



(43) 



Cp = (C I=1 - C I=0 ) /4, 

with C I=1 = -5.021 x 1(T 5 MeV" 2 and C I=0 = -5.714 x 1(T 5 MeV" 2 . 

For the LO2 transfer matrix we use the improved lattice transfer matrix defined in 



(44) 



Mi 



lo 2 



: : exp ^ -H bee a t - ^ ^ /(? 



+ 



8/X 



(45) 



Si,S2,I fti,n2 

where the momentum-dependent coefficient function f(q 2 ) is defined as 

f(q 2 ) = U 1 ex P ~ b Yl ^ " cos ®) 

I 

and the normalization factor /o is determined by the condition 



(46) 



-6^(1 -cos ft) . (47) 

~ q L i 

As in [3| we use the value 6 = 0.6. This gives approximately the correct average effective 
range for the two S'-wave channels when C and Cp are properly tuned. For C and Cp we 
use C I=1 = -3.414 x 10" 5 MeV~ 2 and C I=0 = -4.780 x 10" 5 MeV~ 2 and the relations in 
Eq. (1131 and (l4"4l . More details on these leading-order actions can be found in jl^ . 



B. Next-to-leading-order actions NLOi and NLO2 

For the next-to-leading-order transfer matrices Mnloi an d ^nlo 2 we &dd the following 
local interactions to the leading-order transfer matrices M^o 1 and Mlo 2 - We first start 
with corrections to the leading-order contact interactions. These can be written as 

AV = ±AC:J2p a '' a (n)p a '' a (n):, (48) 

n 

AVp = 1 -ACp :Y,pf' a (n)pf' a (n) : . (49) 

n,I 

At next-to-leading order there are seven independent contact interactions with two deriva- 
tives. These can be written as 

V? = -\c# : ^/> a (n)V 2 /' a (n) :, (50) 



n,l 
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V Ps 2 = - l -Cpj : J^pf ' a (n)vfpf ' a (n) :, 

n,I,l 

^s 2 ,? 2 = -\ c s\f ■ ^2p a s' a {n)V 2 lP i' a {n) :, 

it,S,l 
n,S,I,l 

v {q .s ? = \c [q . s? :j2J2 A sp a s ,a ^J2 A s'pi' a ^) : > 

n S S> 
n,I S S' 

V (iq ,sy k = ~C {iqxS y k : £ £ W [nf' a (n)A,p a ;- a (n) + </(n)A,p at ' a (n) 



(51) 
(52) 
(53) 
(54) 
(55) 
(56) 



The subscripts indicate the continuum limit of the interactions. A detailed discussion of 
the continuum limit for each of these NLO lattice interactions is given in the appendix. 
The V( iqx syk term corresponds with the continuum interaction 



— # 

C( iqxS yk (iq x (ai + a 2 )) ■ k. 



(57) 



We note that this interaction vanishes unless the total spin 5=1. We note also that 

— * 

the continuum limit of the interaction is antisymmetric under the exchange of q and k. 
Therefore the continuum interaction is nonzero only for odd parity channels. Unfortunately 
the lattice interaction V^ qx s)-k does not share this exact t-u channel antisymmetry at nonzero 
lattice spacing. Therefore V^ qx sy k produces small lattice artifacts for S = 1 in even parity 
channels. Fortunately in this case there is a simple way to remove them. We include an 
explicit projection onto total isospin 1=1, 



i=i 

(iqxS)-k 



-jpliqlsyk I 4 : £l < s > 1 

n,l,S,V 



+ 



£i,s,i' 



Iif^{n)A vP f a {n) + </(n)A,p^(n) 
uf/(n)A lip if(n) + uf£(n)A vP f * 



(58) 



n,l,S,l',I 

In the continuum limit the interaction is already pure 1 = 1 due to total antisymmetry and 
so the isospin triplet projection has no effect. At nonzero lattice spacing this projection 
completely eliminates lattice artifacts in the S = 1 even parity channels. 
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FIG. 1: Spherical wall imposed in the center-of-mass frame. 



V. TWO-NUCLEON SCATTERING ON THE LATTICE 
A. Spherical wall method 

We measure phase shifts by imposing a hard spherical wall boundary on the relative 
separation between the two nucleons at some chosen radius -R wa ii [13] . The spherical wall 
removes copies of the interactions due to the periodic boundaries of the lattice. Viewed 
in the center-of-mass frame we solve the Schrodinger equation for spherical standing waves 
which vanish at r = -R wa n as indicated in Fig. [TJ 

In the spin singlet case for values of r beyond the range of the interaction, the spherical 
standing wave can be decomposed as a superposition of products of spherical harmonics and 
spherical Bessel functions. Explicitly we have 

[cos5 L ■ j L (kr) - sm8 L ■ y L (kr)} Y L>Lz (9, 0), (59) 

where the center-of-mass energy of the spherical wave is 

P P 

^ 2 ( 1 

2m m ' 

and the phase shift for partial wave L is Sl- Therefore we know k from the energy E, and 
the phase shift 5l is determined by setting the wavefunction in Eq. ( |59l) equal to zero at the 
wall boundary 

cos5 L ■ j L {kR wa \\) = sin5 L • y L (kR WSb n), (61) 
jL(kR wa a) 



5l = tan 1 



VLikRwaXl) 



(62) 
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On the lattice there is some ambiguity on the precise value of -R wa ii since the components 
of r must be integer multiples of the lattice spacing. We resolve this ambiguity by fine- 
tuning the value of -R wa ii for each standing wave so that 5l equals zero when the particles 
are non-interacting. 

For the spin triplet case, however, spin-orbit coupling produces mixing between partial 
waves L = J — 1 and L = J + 1. If we use the two-component notation, 



Rj-i(r) 
Rj+i(r) 



(63) 



for the radial part of the wavefunction, then we find two sets of standing wave solutions of 
the form 



(X 



kh 



at energy E 1 = (k 1 ) 2 jm and 



m 11 



oc 



k H 7 



AS_ lS in (fcV-^TT + A^J 
A5 +1 sin(^r-^7r + A5 +1 ) 

A 1 /^ sin (k n r - ^tt + Ay_ x ) 
A'f +1 sin (k"r - ^vr + A'/ +1 ) 



(64) 



(65) 



at E n = ) /m. These can be used to derive the phase shifts and and mixing 



angle ej 



tan (— Aj_ x + tan (— Aj +1 + 5j + i) = tan 2 ej, 

tan (-Ay_! + tan (-A^ + 5 J+1 ) = tan 2 ej, 

sin (-A5 +1 + 



A^iitanej = -A 



cos (— Aj_ x + 
sin (-A 7 / +1 + S J+1 ) 



(66) 
(67) 

(68) 
(69) 



■ J+1 COS (-A2.! +*,_!)• 

The phase shifts and mixing angle in Eq. (|66|) and (1681) are at momentum k 1 while the 
phase shifts and mixing angle in Eq. fIBTI) and (I6"9"j) are at momentum /c 7/ . We therefore 
consider only close pairs of values k 1 m k 11 in solving Eq. (I66l) - fl69l) . This can be done 
for example by considering the (n + l) st -radial excitation of L = J — 1 together with the 
n th -radial excitation of L — J + 1. In this scheme we use 



tan (-A5_! + Sj-^k 1 )) tan (-A^ +1 + dj+^k 1 )) = tan 2 [ej^k 1 )] , 
tan (-Aji x + Sj^k 1 )) tan (-A^ +1 + <W^)) ~ tan 2 [e J (k I )\ , 
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(70) 
(71) 



A'^n M*')] = -^ r { r^l' J+ f2 . (72) 

cos (-A J J _ 1 + 

for the phase shifts and mixing angle at k = k 1 , and 

tan (— Aj_ : + 5j-i(k H )) tan (-A 7 J+1 + <$ J+1 (*")) « tan 2 [ej(A; JI )] , (73) 
tan (-Aj 7 _ 1 + 5j_i(A; /J )) tan (-A'/ +1 + 6 J+1 (k n )) = tan 2 [e J (fc // )] , (74) 

cos (-Ay„ 1 + 5j_i(A; JJ )) 
for the phase shifts and mixing angle at k — k 11 . 

For momentum less than or equal to the pion mass all of the mixing angles Ej are 
numerically small. It is therefore convenient to expand in powers of the mixing angle, 

<W^) = A 7 ,-i + : ( A f j{k '\ tut, + 0(e% (76) 
tan(-A J J+1 + 5 J+1 (k 1 )) 



at k = k 1 and 



eAk 1 ) = -^±i sin ( A y +i _ A / + J + 0(4), (77) 



<W*") = A'/ +1 + - , J j{kH ] — , + 0(4), (78) 



tan (-A?..! + <$.,_!(*")) 
5 + i 



ej (*") = -^i sin (A^ - AJ_J + 0(4), (79) 



at fc = k 11 . 



B. LO energy levels and NLO corrections 

In Fig. [2] we show energy levels for spin S = and isospin 1 = 1 using lattice actions 
LOi and LO2. The spherical wall is at radius i? wa ii = 10 + e lattice units where e is a small 
positive number. We use this e notation to make explicit that \ f\ = 10 lattice units is inside 
the spherical wall but all lattice sites with |r| > 10 lattice units lie outside. The solid lines 



indicate the exact energy levels which reproduce data from the partial wave analysis of [40 ]. 
The energy levels for the standard action LOi are about 10% to 15% too low for the 1 So 
states, while the improved action LO2 is correct to within a couple of percent for all 1 S'o 
states. Deviations for the higher partial waves are smaller than 1% for both LOi and LO2. 
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LO 1; 5 = 0,/= 1 



LO 2 ,S = 0,/= 1 



> 




■1 G A 



■ID, 



J z mod 4 



FIG. 2: Energy levels for S = 0, I = 1 using lattice actions LOi and LO2 and a spherical wall at 
radius i? wa ii = 10 + e lattice units. The solid line indicates the exact energy levels which reproduce 
data from the partial wave analysis of 401 ] . 



The energy levels for spin 5 = 0, isospin 1 = 0, and R wa a = 10 + e lattice units are shown 
in Fig. [3] In this case LOi is better for the 1 P\ states and is within 1% of the exact values. 
The LO2 energy levels are further away, though still within 5% for the l P\ states. The error 
for the 1 i 7 3 partial wave is much smaller than 1% for both LOi and LO2. Results for the 
LOi and LO2 in the spin-triplet channels are similar in terms of their relative errors. 

At next-to-leading order we have nine unknown operator coefficients to fit: AC, AC/2, 
C q 2, Cp q 2 } C S 2 q 2, C S 2 tI 2 q 2, C ( g. s) 2, C / 2 j ( ? . g )2, and C!ri_ k . We fit these nine operator 
coefficients using the eight energy levels listed in Table [TT] for i? wa n = 10 + e lattice units, 
as well as the quadrupole moment of the deuteron Qd- The deuteron quadrupole moment 
is a measure of the low-energy strength of S-D partial wave mixing. We could instead use 
the mixing angle E\ at some low momentum scale, but the quadrupole moment is actually 
an easier observable to compute. In this study we incorporate the nine NLO interactions 
using perturbation theory. 
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FIG. 3: Energy levels for S = 0, I = using lattice actions LOi and LO2 and a spherical wall at 
radius i? wa ii = 10 + e lattice units. The solid line indicates the exact energy levels which reproduce 
data from the partial wave analysis of 401 ] . 



The motivation for using perturbation theory for the NLO interactions is that it is com- 
putationally less expensive. Whether or not the NLO corrections can be treated using 
perturbation theory depends on the observable of interest. For example the standard action 
-1O1 has a clustering instability which results in a strong overbinding of the alpha particle 
1^ | . This instability is severe and probably cannot be fixed by introducing NLO corrections 
perturbatively. The problem is that the S'-wave interactions are too attractive at momenta 
q ~ m n , and the S'-wave effective range corrections provided at NLO must be included 
nonperturbatively. 

On the other hand the improved action LO2 resolves this problem using Gaussian smear- 
ing. The Gaussian-smeared interaction is contained in the term V® ^ denned in Eq. (fT51) . 
Since V 2 i s iterated nonperturbatively along with rest of the leading-order action, the 
clustering instability is removed. We use this approach as a general strategy. The first 
attempt is to try a completely perturbative approach for all NLO corrections. If this fails 
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TABLE II: Results for LOi and LO2 and target values 



Spherical wave 


Free nucleons 


LOi 


L0 2 


PWA93 


l l So (MeV) 


0.928 


0.368 


0.418 


0.407 


3 1 S (MeV) 


8.535 


5.838 


6.843 


6.815 


l 3 S(D) l (MeV) 


0.928 


-2.225 


-2.225 


-2.225 


3 3 5(L>)i (MeV) 


8.535 


4.878 


5.430 


5.675 


2 1 P 1 (MeV) 


5.691 


5.755 


5.541 


5.782 


2 3 P(F) (MeV) 


5.691 


5.569 


5.396 


5.584 


2 3 P(F)i (MeV) 


5.691 


5.754 


5.652 


5.753 


2 3 P(P) 2 (MeV) 


5.691 


5.684 


5.558 


5.669 


Qd (fm 2 ) 


N/A 


0.143 


0.276 


0.286 



and some of the interactions must be handled nonperturbatively, then we include those 
interactions in the I/Q 2 /^ 2 term in the leading-order action. 

For each of the nine observables in Table [TT] we compute the derivative with respect to 
each of the nine NLO coefficient operators. By inverting this 9x9 Jacobian matrix we 
find the required values for the operator coefficients needed to match each of the observables 
using first-order perturbation theory. The results for the operator coefficients are shown in 
Table IIHI Although most of coefficients have generally the same order of magnitude and 
sign for NLOi and NLO2, we see that the coefficients in some cases are quite different. 



VI. RESULTS 



We compute lattice phase shifts and mixing angles using spherical walls with radii R wa ii = 
10 + e, 9 + e, and 8 + e lattice units. In order of increasing momentum, the lattice data 
corresponds with the first radial excitation for -R wa n = 10 + e, 9 + e, and 8 + e; second radial 
excitation of -R wa ii = 10 + e, 9 + e, and 8 + e; and so on. The S'-wave phase shifts for LOi 
and NLOi versus center-of-mass momentum pcm are shown in Fig. HI We compare these 



with the S'-wave phase shifts for LO2 and NLO2 in Fig. [51 
are both in good agreement with partial wave results from 



'he NLOi and NL0 2 results 



40j. Systematic errors are 



first noticeable at momenta greater than about 80 MeV and are larger for NLOi. In both 
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TABLE III: Results for NLO operator coefficients 



Coefficient 


NLOi 


NL0 2 


AC (MeV~ 2 ) 


-1.43 x 10- 4 


-1.10 x 10- 5 


ACp (MeV~ 2 ) 


3.42 x 10- 5 


1.03 x 10- 5 


C 2 (MeV~ 4 ) 


1.27 x 10- 9 


-1.01 x 10- 9 


Cprf (MeV" 4 ) 


-5.83 x 10- 10 


-4.07 x 10- 10 


C S 2 :q2 (MeV" 4 ) 


-4.18 x 10- 12 


-1.72 x 10- 10 


C S 2j2 tq 2 (MeV^ 4 ) 


-4.31 x 1(T 10 


-2.81 x 10- 10 


C {q . s) 2 (MeV" 4 ) 


-2.52 x KT 11 


-3.23 x 10- 10 


Cj2 Xq . s) 2 (MeV" 4 ) 


6.24 x 10- 11 


1.43 x 10- 10 


C^syk ( MeV " 4 ) 


1.60 x 1(T 10 


9.81 x 1Q- 11 




FIG. 4: S'-wave phase shifts versus center-of-mass momentum for LOi and NLOi. 



cases the deviations are consistent with higher-order effects such as two-nucleon contact 
interactions with four derivatives. 

The P-wave phase shifts are presented in Fig. [6] and [71 The phase shifts are already not 
bad for LOi and are quite accurate for NLOi. This indicates that at low momenta only a 
small correction is needed on top of P-wave interactions produced by one-pion exchange. In 
the case of LO2 we see that the effect of Gaussian smearing, while useful for S'-wave phase 
shifts, produces an unphysical attraction in each P-wave channel that must be cancelled by 
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FIG. 5: S'-wave phase shifts versus center-of-mass momentum for LO2 and NLO2. 



the NLO2 corrections. For the NLO2 results the residual deviations appear consistent with 
effects produced by higher-order terms such as four-derivative contact interactions. For 1 P\ 
and 3 Po there is some indication of higher- derivative effects for pcm near 110 MeV. 

The .D-wave phase shifts are shown in Fig. [8] and [9j None of the .D-wave data was used 
in the fitting of operator coefficients. For both NLOi and NLO2 results the errors appear 
consistent with effects from higher-order interactions. The NLOi deviations are somewhat 
smaller, though the NLOi and NLO2 deviations appear similar in character. Overall the 
differences among LOi, LO2, NLOi, and NLO2 results for the D waves are smaller than the 
corresponding differences for the S and P waves. This observation is consistent with the 
dominance of the one-pion exchange potential and validity of the Born approximation in 
higher partial waves. 

The mixing parameter e\ in the Stapp parameterization 41[ is shown in Fig. [101 Results 
for LOi and NLOi are on the left, and results for LO2 and NLO2 are on the right. The 
pairs of points connected by dotted lines indicate pairs of solutions at k = k 1 and k = k 11 
for the coupled 3 Si- 3 Di channels. For LOi it is interesting to note that E\ has the wrong 
sign. For both NLOi and NLO2 results the remaining deviations appear consistent with 
effects produced by higher-order interactions. 
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FIG. 6: P-wave phase shifts versus center-of-mass momentum for LOi and NLOi. 



VII. SUMMARY AND DISCUSSION 



We have studied nucleon-nucleon scattering on the lattice at next-to-leading order in 
chiral effective field theory at momenta less than or equal to the pion mass. Throughout 
our analysis we tested model independence at fixed lattice spacing by repeating calculations 
using two different lattice actions. The first leading-order action LOi included instantaneous 
one-pion exchange and same-site contact interactions. The second leading-order action 
LO2 included instantaneous one-pion exchange and Gaussian-smeared interactions. We 
computed next-to-leading-order corrections for these actions and in each case found results 
accurate up to corrections at higher order. In the second paper of this series we use the 
LO2 and NLO2 actions to compute the ground state of dilute neutron matter using Monte 
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FIG. 7: P-wave phase shifts versus center-of-mass momentum for LO2 and NLC>2- 

Carlo. This is done using the auxiliary-field transfer matrix method introduced in fl3 |. 

Overall we find that the Gaussian-smeared actions LO2 and NLO2 are more accurate 
than the standard actions LOi and NLOi. This can be seen most easily from the size of 
the required NLO corrections in Table [Til For P-wave interactions, however, we see from 
comparing Fig. [6]and Fig. [7]that the standard actions LOi and NLOi are more accurate. In 
future studies it would be useful to try to find an improved LO action with accurate S"-wave 
phase shifts and smaller P-wave attraction without inducing sign problems in Monte Carlo 
simulations. 

The results presented here can be extended to higher orders in chiral effective field theory. 
If we continue with the low cutoff momentum A w then the NNLO two-pion exchange 

potential can be expanded in powers of q 2 / (4m^) in the same manner as the NLO two-pion 
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FIG. 8: D-wave phase shifts versus center-of-mass momentum for LOi and NLOi. 



exchange potential. This expansion yields local operators which are either renormalizations 
of local operators with zero or two derivatives, local operators with four derivatives, or local 
operators with more than four derivatives. The operators with four derivatives should be 
treated in the same manner as local four-derivative operators appearing at N 3 LO in the 
usual chiral power counting. 

In our analysis we have ignored small effects that appear in extraneous channels due to 
broken rotational invariance on the lattice. For example the mixed 3 5'i- 3 D 1 spherical waves 
each have some small admixture of an unphysical 3 D 3 component. In our analysis of E\ we 
measured the ratio of 3 S\ and 3 D\ components and neglected this small 3 D% component. 
In the future we might prefer a more ambitious approach which explicitly removes the 3 D 3 
component. This requires including local interactions which are invariant under the cubic 
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FIG. 9: .D-wave phase shifts versus center-of-mass momentum for LO2 and NLO2. 



rotational group SO(3, Z) but not invariant under the full SO (3) symmetry. For example 
in addition to the SO(3)-invariant interaction 

\- E E A S P f a (fi) £ ^s> P i' a (n) :, (80) 

n S S' 

we should also include a small contribution from the SO(3, Z)-invariant but SO(3)- 
noninvariant operator 

1 -:J2J2 A sPs' a (n)A sP ^ a (n):. (81) 

n S 
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FIG. 10: E\ mixing angle for LOi and NLOi on the left, LO2 and NLO2 on the right. 
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APPENDIX A: CONTINUUM LIMIT OF THE NLO INTERACTIONS 
1. One-nucleon matrix elements 

The matrix elements of the local densities for single nucleon states are 

(p,i,i|/'°(n) \p',i',j') = dwS^'-^, (Al) 

{p,i,j\pf' a (ft) \p',i',j') = 6* [Tt] jjt eW-*)*, (A2) 
(P,hj\p a s' a (n) \p',i'J') = [0s]u,6 jj .e i V-*)* (A3) 
(P^j\p a sf(n) \p',i',f) = [a s ] zz , [n\ ir (A4) 
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The matrix elements for derivatives of these densities are found by taking f(n) = e 1 ^' 
and computing 

A,/(n) = j®'-*)-* x ^ (-ly^eW-fr*, (A5) 

V 2 J(n) = e^'-fr* x p(p'-f> f + e -i(F-f)-i _ 2] . (A6) 
In the continuum limit \p'\ and \p\ in units of inverse lattice spacing are small and so 

\ E (-lY^e^'-^^iip't-Pi), (A7) 
^1,^2,^3=0,1 

e i(p-'- P -)-i + e -i(p-'- P -)-J _ 2 _> _ (p| _ . (A8) 
The matrix element for the SU(4)-invariant current density is 

(p,i,j\ nf ' a (n) |p',i',/> = Su^e^'-^ x 1 £ ( _i)^ + i e ^'-V^^'). (A9) 

^1,^2,^3=0,1 

Similarly the matrix element for the spin current density is 

(p,i, j\Tlf/{n) \p',i',f) = [vs^foeW-fi-* x \ £ (_i)^+i e ^-V^-D. ( A 10) 

^1,^2,^3=0,1 

The other current densities are not needed for the NLO interactions. In the continuum 
limit 

I ^ (_ 1 )^+i e ^'- e -^--(-0 _ j (pj + Pl) . ( A11 ) 

^1,1^2,^3=0,1 

2. Two-nucleon matrix elements 

The tree-level amplitude for two-nucleon scattering consists of contributions from direct 
and exchange diagrams. However for bookkeeping purposes we label the contact interactions 
according to the tree-level amplitude for scattering of two distinguishable nucleons. We 
imagine that one nucleon is of type A and the other nucleon is of type B. The interactions 
include densities and current densities for both A and B. For example the SU(4)-invariant 
density and current density become 

p a ^ a {n) -> p a A' aA (n) + p a *' aB (n), (A12) 

Uf A (n) -> uf A ' aA (n) + Tlf B ' aB (n). (A13) 
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Let the incoming momenta be p^ and pf and the outgoing momenta be pf and pf . The 
t-channel momentum transfer is 



-> ->A -^B , ->B 

Q = Pi ~Pf = ~Pi +Pf, 



and the w-channel exchanged momentum transfer is 

k=P l A -pf = -p l B +pf- 
For these incoming and outgoing states the amplitudes for AV^ are 

A (AF) = AC, 

^(Ay 2 ) = AC /2 E r f r f- 
i 

In the continuum limit the amplitudes for are 

A(V q2 )^C q2 q 2 , 
A (y 2)(?2 ) -> C^g 2 E r/rf, 

^(y^^c^^E^f, 

s 

A (V S 2j2tf) -> C S 2j2^q 2 ^ °i°~f E T ^ T ^' 

-4 (y^.s) 2 ) -> C( 9 -5) 2 X] 5^ 9s ,£r f'> 

5 5' 

A (Vp^.sy) -> C/2 )(9 . 5) 2 E r/rf E g s <7g E 9s'<rf , 
-4 (V(i 9 xS)-fc) — » iC(j qxS )-k eisi'Qi (v A + v B ) s . 

l,S,l' 



(A14) 



(A15) 



(A16) 
(A17) 



(A18) 
(A19) 

(A20) 

(A21) 

(A22) 

(A23) 

(A24) 



Most of these are straightforward, however the spin-orbit amplitude A (V(iqxs)-k) requires 
some derivation. In terms of the incoming and outgoing momenta, 

A (V(iqxS)-fc) 

- E ^ \(P? + Pf), o» 8 (pf - p B ) v + (pf + pf ), a A s (pf - p}) p 



+ 2<w* E W + p})i °a (pf - p?) P + (pf + pf W W - 



(A25) 
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We note that 



E (Pf + Pfh-s (pf ~ P? ) v = ~ E \(pf + pf)rf (pf - pj) v 



i,s,v 



J2^\(2pt) l ^(pt-Pf) l ]- (A26) 



Similarly 



E £ « (pf + pf W (pf - pf), = - E [(*?), (pf - pf),. 

Z,S,J' Z,S,Z' 

E ««■ (pf + p; 1 ), (pf - pf ), = - E ««■ W W (pf - pf),-' 



Z,S,Z' 



Z,S,Z' 



E ««- ( P f + p?W (pf -?/)„ = -E £ « [(*?), S (Pf - Pf ), 

l,S,l' l,S,l' 

The sum of the four terms in Eq. flA26j) - f)A29l) gives 



E £isi> 

l,S,V 

which equals 



We therefore find 



E eisv [-2ha§ qy - 2kia§q v ] . 



i,s,v 



A (V( iqxS yk) = -iC {ig>< syk E £lsl ' [ k i a sQi' + hvsQi'] 

i,s,v 



(A27) 

(A28) 
(A29) 

(A30) 
(A31) 



(A32) 
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